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ABSTRACT 


During  this  program,  the  specifications  for  the  ANDVT  HF  modem 
have  been  refined  and  detailed  evaluation  and  simulation  of  the  new 
technical  features  within  this  modem  have  been  conducted.  These 
include  a multiple -tone  signal  detection  format  with  an  adaptive 
threshold,  a multiple-tone/multiple-stage  Doppler  estimation  algo- 
rithm, a matched  filter  frame  estimation  algorithm  utilizing  PN 
correlation  properties,  a low-rate  error-correction  coding  approach 
for  protection  of  the  KG  ' ync  sequences,  an  error-correction  coding 
approach  specifically  4esigned  to  protect  the  critical  speech  para- 
meters, use  of  soft-decision  (channel  measurement)  information  ob- 
tained from  the  demodulator,  and  decision-directed  Doppler  tracking 
utilizing  information  from  all  data  tones . The  analytical  and 
simulation  results  provide  the  desirable  result  that  the  preamble 
can  be  successfully  received  at  a lower  SNR  than  is  required  for  the 
reception  of  high-quality  2400-b/s  digitized  voice.  On  a previous 
effort,  high-quality  digitized  voice  reception  was  obtained  over  a 
time-varying  multipath  HF  fading  channel  model  at  a signal-to-noise 
ratio  per  information  bit,  E^/Nq,  in  the  range  of  11  to  14  dB.  The 
preamble  functions  evaluated  during  this  program  typically  required 
several  dB  lower  signal-to-noise  ratio  than  is  required  for  success- 
ful reception  of  digitized  voice. 
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SECTION  1 


INTRODUCTION  AND  RECOMMENDATIONS 


This  report  represents  a continuation  of  the  previous  work 
done  by  CNR  [1.1]  and  the  Naval  Research  Laboratory  [1.2]  on  the 
design  and  evaluation  of  an  HF  modem  for  the  advanced  narrowband 
digital  voice  terminal  (ANDVT) . The  emphasis  in  this  report  has 
been  on  computer  simulations,  whereas  the  previous  report  [1.1] 
emphasized  the  specifications  and  theoretical  analysis  of  the 
signaling  formats.  As  a result  of  this  effort,  a high  degree  of 
confidence  can  be  placed  on  the  advanced  modem  techniques  speci- 
fied. Table  1-1  summarizes  the  technical  advances  in  HF  modem 
design  which  have  resulted  from  this  and  the  previous  study  [1.1] 
sponsored  by  the  Navy. 


The  recommended  format  is  illustrated  on  Figure  1.1. 

Signal  detection  and  Doppler  estimation  is  accomplished  in 

approximately  0.355  sec  by  first  detecting  the  presence  of  a 

signal  and  then  estimating  the  Doppler  offset  which  may  b:e  as 

1-2 

large  as  + 75  Hz.  A probability  of  detection  of  below  10  and 
negligible  probability  of  a false  alarm  is  achievable  by  this 
format.  The  Doppler  estimation  error  is  well  below  0.5  Hz. 
Reference  [1,1]  and  the  simulations  presented  in  Section  3 of 
this  report  are  in  excellent  agreement  and  indicate  that  the 
above  performance  results  can  be  achieved  at  a signal-to-noise 
density  ratio*  of  AO  dB-Hz,  which  corresponds  to  a signal-to- 
noise  ratio  measured  in  a 2400-Hz  bandwidth  of  only  6.2  dB. 


The  signal  to  noise  density  ratio  P/Nq  is  related  to  the  E^/Nq 
for  a 2400-b/s  voice  modem  by: 


ft) 


voice 


+ 33.8  dB 


1-1 


1 U Bl 

o 

8 § 3 

o 5 § 

g 5 § 

u w w 

M Q O 

02  I 

g t s 

O H 

3 w s 

*J  W Q 


h O Q 
M H W 


0 3 

u < 


rt  o 


Pd  H 

o H 


DC  Dd 

w o 


ALL  DATA  TONES 


Figure  1.1  A Recommended  Half-Duplex  Format 


noise  ratio  measured  in  a 2400-Hz  bandwidth  of  only  6.2  dB. 


Following  the  signal  detection  and  Doppler  estimation  format, 
a frame  estimation  signaling  format  composed  of  wideband  («  1440Hz) 
pulses  forming  three  PN  sequences,  each  of  length  31  pulses,  is 
transmitted.  This  framing  technique  is  shown  in  Section  2 of  this 
report  to  be  superior  to  the  Carley  framing  technique  specified  in 
[1.1J,  The  simulation  of  this  framing  format,  documented  in  Sec- 
tion 3,  indicates  that  framing  errors  of  below  1 ms  are  achieved  at 
signal-to-noise  density  ratios  of  40  dB-Hz  for  realistic  2-  and 
6-path  HF  channel  models.  It  should  be  noted  that  the  simulations 
reported  on  in  Section  3 include  the  entire  sequence  of  events,  i.e., 
signal  detection,  Doppler  estimation,  and  frame  estimation,  in  a 
manner  suitable  for  hardware  implementation. 

Following  the  PN  frame  estimation  format,  the  39-tone  four- 
phase  PSK  modem  begins  transmission.  The  next  0.6525  second  is 
devoted  to  the  reception  of  the  KG  synchronization  sequence.  Ten 
Golay  codes,  with  eighth-order  diversity,  are  transmitted  over 
25  frames.  The  first  22.5-ms  frame  is  used  as  a phase  reference 
and  the  last  three  frames  are  used  for  soft  decoding  of  these 
Golay  codes.  Since  the  actual  data  rate  is  very  low,  i.e., 

120  bits  are  transmitted  in  0.6525  sec,  results  are  obtained  in 
terms  of  the  signal-to-noise  per  tone  E^/Nq,  which  is  closely 
related  to  the  signal-to-noise  ratio  measured  in  a 2400-Hz  band- 
width, i.e.,  | . For  the  two-path  HF  channel  modem,  the 

\ vAfoice 


*As  noted  in  Appendix  B,  the  signal-to-noise  ratio  per 
is  related  to  P/Nq  by: 


+ 33.4  dB 


tone  Et/Nq 
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simulation  results  in  Section  4 indicate  that  a probability  of 

_3 

a block  error  on  120  bits  of  below  10  is  achievable  at  a 

p 

E.f/No  ■ 8 dB  or  — = 41.4  dB-Hz.  For  the  six-path  split  ray  model, 
a probability  of  a block  error  of  below  10  ^ is  achievable  at 

p 

E^/N  = 3 dB  or  —*36.4  dB-Hz.  These  results  are  consistent  with 
TO  N0 

our  design  goals  of  requiring  the  preamble,  including  the  KG  sync 
sequences,  to  operate  at  an  SNR  below  that  required  for  acceptable 
voice  c ality. 

Tapes  of  voice  transmission,  based  on  the  NRL  LPC-10,  which 

were  made  during  the  previous  effort  [1.1]  indicated  that  an  E^/Nq 

of  14  dB,  or  —-*47.8  dB-Hz,  is  required  for  the  two-path  model, 

, nQ  P 

and  an  E,  /Nn  of  11  dB,  or  — = 44.8  dB-Hz,  is  required  for  the 
D U Nq 

six-path  split  ray  model.  Thus,  the  entire  preamble,  signal  pre- 
sence, Doppler  estimation,  frame  estimation,  and  KG  synchronization 
can  be  received  at  several  dB  lower  than  that  required  for  intelli- 
gible voice  receptions.  Unfortunately,  no  new  voice  tapes  have 
been  made  during  this  effort  since  the  FORTRAN  version  of  the 
latest  LPC-10  format  supplied  by  the  Government  was  inoperative. 

During  the  reception  of  the  KG  synchronization  sequence  and 
the  reception  of  voice  transmission,  Doppler  and  frame  tracking 
is  required.  Simulated  results  for  the  new  Doppler  tracking  algo- 
rithm are  given  in  Section  3 of  this  report  and  compared  to  the 
more  conventional  single-tone  tracking  algorithm.  The  new  algo- 
rithm, which  is  based  on  using  all  of  the  modulated  data  tones, 
offers  better  performance  than  the  single-tone  Doppler  tracking 
loops  and  has  the  added  advantage  that  no  additional  bandwidth  or 
power  is  required  for  Doppler  tracking.  Frame  tracking  is  achieved 
by  using  early  and  late  integration  gates  within  an  empty  tone  slot 
"slot  sync",  which,  incidentally,  is  the  only  conventional  HF  modem 
technique  used  in  this  modem. 
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SYNCHRONIZATION  PREAMBLE  PERFORMANCE 
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2.1  Introduction 

In  this  section  we  make  a careful  examination  of  the  synchro- 
nization preamble  performance.  Recall  from  earlier  work  that 
the  correlation  peaks  at  the  matched  filter  outputs  and  other 
properties  of  sync  preamble  performance  (e.g.,  effect  of  Doppler 
offset)  are  all  derived  from  consideration  of  a transmitted  sync 
preamble  of  the  form 

M-l 

zn(t)  « X V,  P(t  - iA)  (2.1) 

U i=0  1 

where  the  discrete  elements  [y^,  i = l,  ...»  M-l}  constitute  the 
synchronization  sequence  and  are  chosen  for  their  correlation 
properties.  Because  there  is  no  periodicity  in  (2.1)  the  discrete 
matched  filter  response  to  the  synchronization  sequence  in  (2.1) 
has  a special  form  and  is  referred  to  as  the  aperiodic  correla- 
tion function. 

For  completeness  and  future  comparisons,  the  aperiodic  cor- 
relation function  is  developed  in  Appendix  A [see  Eq.  (9)  for 
explicit  representations] . In  this  section  we  approach  the  problem 
from  a more  general  point  of  view  which  incorporates  features  of 
both  the  periodic  and  aperiodic  correlation  functions.  This 
topic  is  treated  in  Section  2.2  which  has  as  its  main  objective  an 
ambiguity  function  analysis  of  the  performance  of  a sync  preamble 
using  PN  sequences  and  an  analysis  of  the  edge  effects  due  to 
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non-periodicities  at  the  beginning  and  end  of  the  preamble.  In 
both  cases  we  are  interested  in  the  magnitude  of  sidelobes  in  the 
matched  filter  output;  in  particular,  we  are  interested  in  the 
extent  to  which  the  sidelobe  levels  are  modified  in  the  presence 
of  a channel  Doppler  shift. 

In  Section  2.3  we  address  the  problem  of  evaluating  the  rms 
estimation  error  associated  with  the  determination  of  frame  synchro- 
nization. The  signal  processing  is  the  same  as  that  described  in 
an  earlier  report  [2.1] at  least  up  to  the  output  of  the  discrete 
matched  filter.  Samples  of  a PN  sequence  correlation  function, 
distorted  by  channel  and  noise,  are  available  at  that  point.  To 
this  collection  of  samples  we  now  apply  a selection  rule  whereby 
one  selects  the  maximum  energy  sample  in  a segment  of  the  discrete 
matched  filter  output  that  is  equal  in  duration  to  one  period  of 
the  sync  preamble.  In  Section  2.3  we  analyze  this  situation  and 
present  results  which  demonstrate,  in  terms  of  rms  estimation  error, 
that  excellent  performance  is  obtained. 

2 . 2 Ambiguity  Function  Analysis 

In  this  section  we  perform  a rather  general  ambiguity  func- 
tion analysis  of  the  sync  preamble.  We  recall  that  we  are  inter- 
ested in  the  sidelobe  level  in  the  output  of  the  matched  filter 
receiver.  With  no  Doppler  shift  present,  the  matched  filter  output 
corresponds  to  the  correlation  function  of  the  discrete  sequence 
used  to  generate  the  sync  preamble.  When  there  is  a Doppler  off- 
set, the  output  corresponds  to  the  discrete  ambiguity  function  of 
the  preamble  and  degradation  can  result.  It  is  one  purpose  of 
this  section  to  determine  that  degradation. 
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Generally,  each  of  the  functions  above  has  a different 
value  depending  on  whether  the  correlator  acts  to  calculate  the 
periodic  or  aperiodic  correlation  function.  When  PN  sequences 
are  used  one  must  be  concerned  with  the  periodic  correlation 
function.  However,  transmission  of  a finite  duration  preamble 
necessitates  examination  of  edge  effects  at  the  beginning  and 
end  of  the  preamble.  This  topic  is  also  treated  in  this  section. 

We  can  generally  refer  to  the  sync  preamble  as  repeated 
versions  of  a fundamental  preamble  element 

M-l 

zQ(t)  = 1 yjPC*  “ *-£)  (2«2> 

i=0 

If  we  make  the  following  definitions: 

P s length  of  period  (integer  multiples  of  A) 

N s number  of  repetitions 

we  can  write  the  total  transmitted  sync  preamble 

N-l 

z(t)  - V zn(t  - j PA)  (2.3) 

j-0 

Thus,  we  see  that  the  total  sync  preamble  consists  of  repeated 
versions  of  z^Ct)  with  zeroes  between  (P  > M)  or  directly 
abutting  (P  * M) . It  is  convenient  to  use  the  representation 

P - qM  + n ^ (2.4) 

where  q and  n are  integers  with  the  restrictions  q > 1 and 
0 < n < N.  The  fact  that  q can  not  be  equal  to  zero  rules  out  the 
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possibility  of  overlap  which  would  be  equivalent  to  a redefini- 
tion of  the  sequence. 

It  is  helpful  to  write  out  (2.3)  in  a way  that  exhibits  the 
dependence  on  the  discrete  sequence  i = 0>  M - 1}.  From 

(2.2)  we  have 

M-l  N-l 

z(t)  - I I y,  p[t  - (i  + jP) A]  (2.5) 

i=0  j=0  1 

We  are  interested  in  evaluating  sync  preamble  performance 
in  the  presence  of  a Doppler  shift  of  F Hz.  Thus,  the  received 
signal  at  the  input  of  the  received  pulse  matched  filter  is 
obtained  by  multiplying  (2.5)  by  exp  C j 2wFt} . After  doing  this 
and  rearranging  we  have 

.(«>  - Y "fy.  eJ2"<i+kP>FAp[t-(i+kPUl  (2-6) 
i=0  k=0 

where 

p(t)  = p(t)ej2,,Ft  (2.7) 


Thus,  a comparison  of  (2.6)  with  the  result  that  is  obtained  when 
F = 0 indicates  a multiplication  of  the  original  coefficients  by 
a cissoid  and  the  pulse  modification  in  (2.7). 

It  is  our  goal  to  determine  the  output  of  the  receiving  pulse 
filter  which  has  an  impulse  response  q(-t),  i.e.,  is  ideally  matched 
to  some  pulse  shape  q(t).  Denoting  this  response  w^,(t),  we 
compute 
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W^t)  - z(t)  ® q(-t) 


(2.8) 


to  obtain 

„ (t)  - ^ ^ v.  ej2”(i+kI')fA  r[t  . (i+kp)A,F]  (2.9) 

1 i-0  k=0 

where 

r(t,F)  - p(t)  ® q*(-t) 

=?  J p(r)q*(T  - t)e+*^irFT  dT  (2.10) 

Because  the  Doppler  offset  F is  small  and  causes  a negligible 
shift  in  the  spectrum  of  p(t)  [see  (2.7)],  we  can  write 


Thus,  we  deal  with  the  matched  filter  output  samples 


w - ^ % v.  ej2nF(i+kP)A  r ( [m  - (i+kP)]A}  (2.15; 

m k=0  i=0  1 

On  the  inner  sum  the  only  contribution  occurs  when  i = m - kP, 
and  the  contribution  takes  place  only  if  the  inequality 


0<m-kP<M  - 1 


To  simplify  matters  wti  use  the  function 


(2.16 


[a,b] (n)  « 


a c n < b 


otherwise 


(2.17 


and,  accordingly,  reduce  the  double  sum  in  (2.13)  to  a single  sum 


%-L  VkP  eJ2"Fm*[0,  M-lKm-kP) 

k=0 


(2.18 


It  is  a simple  matter  to  determine  the  form  of  w for  the  L'th 
r m 

period;  specifically  we  have 


w = 
m 


ym-(L-l)P 


J27rFmA 


(L-l)P  < m < (L-l) P + M-1 


(L-l) P +M-1  < m < LP 


The  representation  in  (2.18),  however,  is  simpler  to  apply 
analytically. 


The  response  of  an  (FIR)  filter  of  length  M - 1 to  the  dis  - 
crete input  in  (2.18)  is  given  by 


M-l 

i«o 


hjW 


n-i 


(2.20) 


For  the  particular  application  we  use  a discrete  matched  filter 


hi  - <2-21> 

The  difference  between  this  choice  and  that  used  in  Eq.  (4)  of 
Appendix  A is  that  here  we  have  made  the  minimum  shift  to  the 
right  required  to  guarantee  causality.  We  have 


N-l  M-l  ,u 

y"  ‘ Jo  Vi-kryM-l-ie  ( } t°.  M-l]  (n-l-kP)  (2. 

k=0  i=0 


With  the  change  of  variable  q - n-i  we  have 


N-l 

- I . 

n k=G  o=n-(M-l) 


n 

V 


yq-kPY*M-l)-t q-n  ej2”FqA  «-D(q-kP) 

(2.23) 


We  shall  find  considerable  use  for  this  expression. 
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One  reduction  of  interest  occurs  when  we  consider  just  one 
period,  i.e.,  when  we  set  N“l,  we  obtain 


w 


n 

V 


n 


n-(M-l)Vtl  r<M-D+q-n 


J2»FqA  [0>  M.1J(q) 


(N-l) 


(2.24) 


For  convenience  in  comparison  it  facilitates  matters  to  evaluate 
(2.24)  at  n ■ p + (M-l).  We  get 


Wp+(M~1) (F> 


M-l 


V 

q»p 


y y * 

q “q-p 


ej2ffFqA 


(2.25) 


which,  unlike  (2.24)  which  peaks  at  q - M-l,  peaks  at  p*0.  This 
result  agrees  exactly  with  an  earlier  result  [2,1,  p.  2-84],  In 
<2.24),  however,  we  have  the  peak  explicitly  taking  place  at  n»M-l 
which  makes  this  expression  more  suitable  for  our  present  purposes, 
especially  when  it  involves  keeping  track  of  the  timing. 


Using  the  general  expression  (2.23)  we  have  calculated  the 
ambiguity  functions  for  several  cases  of  interest.  The  matched 
filter  output  for  a length-31  PN  sequence  is  shown  in  Fig.  2.1  for 
Doppler  shifts  of  0.0,  2.0,  and  10.0  Hz.  The  leading  portion  of 
each  curve  illustrates  the  forementioned  edge  effect.  The  right- 
hand  section  shows  the  internal  structure  of  the  matched  filter 
output  when  total  periodicity  applies;  in  this  segment  the  peak- 
to-sidelobe  ratio  takes  on  the  value  M (where  M is  the  sequence 
length)  when  the  Doppler  shift  F is  zero. 
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Magnitude  of  Discrete  Matched  Filter  Output  for  Length-31  PK  Sequence  (If  sequence 
is  repeated  n times  the  output  consists  of  segment  A followed  by  n-1  repeats  of 
segment  B followed  by  a time  reverse  of  segment  A.)  p.  and  are  the  peak-to- 
sidelobe  ratios  m the  respective  segments.  a b 


For  a Doppler  shift  of  2 Hz  (the  maximum  that  need  be 
allowed),  we  see  that  p^,  the  si  delobe  level  in  the  region 
of  interest,  has  a value  in  excess  of  21  for  the  31  element 
sequence.  This  degradation  does  not  seriously  degrade  the 
performance  probabilities  relative  to  peak  detection. 

For  purposes  of  comparison  we  have  also  provided  data  for 
the  F ■ 10  Hz  Doppler  shift  in  Fig.  2.1.  Here  we  see  that  the 
sidelobe  ratio  has  degraded  from  31  to  5.61  in  the  region  of 
periodicity  whereas  at  the  edges  it  has  only  degraded  from  6.20 
to  5.28,  i.e.,  the  two  sidelobes  are  of  comparable  level.  We 
re-emphasize,  however,  that  this  Doppler  shift  far  exceeds  the 
shift  after  Doppler  correction. 

For  purposes  of  completeness,  the  edge  effects  for  various 
sequences  of  lengths  3 through  7 are  shown  in  Fig.  2.2.  It  can  be 
seen  that  the  degradation  can  be  quite  substantial  (e.g.,  M ■ 11). 
For  this  reason,  gating  out  of  the  areas  at  the  edges  is  recom- 
mended generally. 

A particularly  interesting  sequence  is  shown  in  Fig.  2.3. 
This  sequence  is  the  only  known  optimal  binary  sequence.  Its 
sidelobes  are  zero  when  there  is  no  Doppler  offset,  and  the  peak- 
to-sidelobe  ratio  is  80  when  F ■ 2.0  Hz.  Despite  its  extremely 
favorable  sidelobe  properties,  this  sequence  is  not  useful  in 
our  current  context  as  we  shall  see  in  Section  2.3. 
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Figure  2.2  Illustration  of  Edge  Effects  for  PN  Sequence 
Correlation  Functions 
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2 . 3 Evaluation  of  Frame  Sync  Estimation  Error 

This  section  addresses  the  problem  of  evaluating  the  estimation 
error  associated  with  the  determination  of  frame  synchronization. 

The  signal  processing  is  exactly  the  same  as  that  described  in  an 
earlier  report  [2.1]  at  least  up  to  the  output  of  the  discrete 
matched  filter.  Recall  that  samples  of  the  PN  sequence  correlation 
function,  distorted  by  channel  and  noise,  are  available  at  that 
point.  To  this  collection  of  samples  we  now  apply  a selection 
rule  whereby  one  selects  the  maximum  energy  sample  in  a segment 
of  the  discrete  matched  filter  output  that  is  equal  in  duration 
to  one  period  of  the  sync  preamble.  This  section  provides  an  analysis 
of  this  situation  and  presents  results  which  demonstrate,  i terms 
of  rms  estimation  error,  that  very  good  performance  is  obtained.  The 
processing  is  indicated  schematically  in  Figure  2.4. 

Our  previous  calculations  implicitly  assumed  ideal  sampling 
at  the  pulse  repetition  rate,  hence  our  results  which  depended 
on  the  discrete  autocorrelation  function  of  the  sequence  (PN) 
that  modulates  the  pulse  train.  In  practice,  the  synchronization 
pulse  duration  will  be  some  multiple  of  the  basic  sampling  interval 
A,  e.g.,  IcqA.  In  such  a case,  because  we  cannot  a priori  separate 
the  peaks  from  the  sidelobes,  we  must  correlate  kA  sets  of  N 
variables  corresponding  to  the  samples  taken  kQA  seconds  apart. 

The  three  sets,  of  course,  arise  from  the  three  possible  "guesses" 
of  the  "ideal"  sampling  times.  The  best  choice  of  sample  point 
sets,  i.e.,  the  one  that  provides  the  sample  closest  to  the  peak, 
is  marked  with  arrows  in  Figure  2.5.  Note  that  the  choice  of  the 
best  set  does  not  provide  a sample  corresponding  precisely  to  the 
peak;  however,  the  sample  corresponding  to  the  peak  will  be  in 
error  by  only  one  sampling  interval  at  most.  In  order  to  simplify 
our  calculations,  it  is  assumed  in  the  remainder  of  this  develop- 
ment that  the  peak  value  can  be  obtained  precisely. 
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Selecting  the  peak  value  according  to  the  rule  described 
above  is  mathematically  equivalent  to  simply  choosing  the 
maximum  of  the  sample  points. 

For  the  purposes  of  our  analysis  we  view  the  time  axis  as 
being  divided  into  two  mutually  exclusive  sets: 

(1)  Received  pulse  (and  multipath)  interval 

(2)  Noise  interval. 

The  first  region  is  occupied  with  (multipath-distorted)  signal 
energy.  The  second  region,  because  of  the  extremely  low  side- 
lobes,  is  taken  to  be  occupied  by  noise  energy  only.  We  ignore 
fluctuation  in  estimated  frame  boundary  due  to  multipath  since 
multipath  spread  is  small  compared  to  frame  duration.  When 
pulses  overlap  due  to  multipath,  the  peak  can  be  reduced  and  in 
the  limit  of  overlap  one  must  consider  a single  path.  Since  we 
calculate  as  a function  of  non-overlapping  pulses,  we  span  the 
calculation  range  of  interest.  For  those  errors  caused  by  noise, 
the  characterization  of  the  pulse  with  a single  value  corresponds 
to  a worst-case  calculation. 

The  notation  for  this  problem  is  shown  in  Figure  2.6. 
Generally,  the  multipath  components  are  assigned  the  subscripts 
|n^;  i~l,  ...,  L|.  Defining  the  multipath  centroid 


L 


the  error  under  consideration  is  taken  to  be  the  weighted  sum 
over  all  N sample  slots 
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(delay) 
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x = nA 
n 


TOTAL  OBSERVATION  INTERVAL  = NA 


Figure  2.6  Nomenclature  and  Notation 
for  3 -Path  Channel 
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(2.27) 


N 


Z (v*)2  pj 

J-i 


where  is  the  probability  that  the  energy  in  the  sample  slot 
at  Xj  exceeds  the  energy  in  all  the  other  slots.  This  sum  can 
be  divided  up  in  the  following  way: 


■ Z (v*)2  pj + ZKr*)2  pi 

j=i  i-i 

tip  n2,...»np 


(2.28) 


Note  that  the  first  sum  in  (2.28)  is  over  those  slots  that  contain 
no  signal  energy;  i.e.,  it  represents  a measure  of  the  error  due 
to  noise  only.  The  second  sum  in  (2.28)  represents  an  error  arising 
from  the  combined  effects  of  noise  and  fading.  Recalling  that 
the  word  "selection"  will  be  used  to  indicate  choice  of  the 
sample  slot  with  maximal  energy,  we  make  the  following  definitions: 


1)  p is  defined  as  the  probability  that  a sample 
slot  containing  noise  energy  only  will  be 
selected .* 

(2)  q is  defined  as  the  probability  that  a sample 
slot  containing  noise  energy  and  signal  energy 
will  be  selected. 


elicit  in  this  development,  of  course,  is  our  assumption  that 
the  noise  samples  are  independent  and  identically  distributed. 

This  assumption  is  not  strictly  true,  of  course,  because  of 
smearing  by  the  pulse  filters,  but  is  in  keeping  with  our 
commitment  not  to  utilize  assumptions  that  favor  system  operation. 
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Using  the  definitions  in  (2.28)  we  have  , after  some  rearranging, 


t My 

j-1 


+ (q-p) 


(2.29) 


In  the  absence  of  noise,  p=0  and  q=l/L  (recall  the  assumption 
of  equal  strength  paths).  Hence  we  obtain 


e 


2 


(2.30) 


as  the  mean-squared  error  in  the  no-noise  case.  In  (2.30) 


(2.31) 


is  a measure  of  multipath  spread  (multipath  spread  = 2am) . 


Returning  to  (2.29)  we  have 
. N 

e2  - P ( xj“*)2  + L ( q-p  ) am  (2.32) 

Some  thought  about  the  second  term  in  this  expression 

2 

indicates  that  it  is  just  slightly  less  than  crm  for  all 
practical  situations.  In  the  sequel  we  refer  to  it  as  the 
multipath  error  and  neglect  it.  It  can  always  be  reasonably 
approximated  by  (2.30)  once  the  multipath  spread  of  the  particular 
channel  under  study  is  known.  The  remaining  error  term 
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(2.33) 


we  refer  to  as  "the  error  in  excess  of  multipath  spread"  and 
calculate  below. 

To  evaluate  the  sum  in  (2.33)  we  recall  that  • jA  and 
assume  for  simplicity  that  x»mA  where  m is  an  integer  satisfying 
the  inequality  1 < m £ N.  This  assumption  at  worst  introduces 

an  en- :*r  less  than  A/2.  Applying  the  identities 

22 1 " f <m  + i) 

i=l 

and 

m 

Z i2  - | (-  + 1X2-  + 1) 

i*=l 

some  manipulation  yields 

e2  -pA2  [£  + 4+|  + *0o]  (2-36) 

where 

f(m)  ■ mN  (m  -N  - 1)  (2.37) 

Equation  (2.37)  is  a function  of  the  random  variable  m which 
represents  the  centroid  location  within  the  observation  interval. 


(2.34) 


(2.35) 
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It  is  entirely  reasonable  to  assume  that  the  signal  centroid, 
and  thus  m,is  uniformly  distributed  throughout  the  observation 
interval . This  yields 

f - - \ N (N+l) (N+2)  (2.38) 

as  the  mean  value  of  f(m).  For  the  average  error  we  obtain 
2 

e2  - N(N2-l)p  (2.39) 

It  remains  to  evaluate  the  probability  p that  a sample 
slot  containing  noise  energy  only  will  be  the  one  selected. 

Such  a slot  must  contain  more  energy  than  all  the  other  slots ; 
this  includes  slots  with  signal  and  noise  energy  as  well  as 
other  slots  with  noise  energy  alone.  In  the  noise  only  slots, 
we  assume  that  z^,  which  represents  the  magnitude  squared  value 
of  the  signal,  is  distributed  according  to  the  probability  density 
function  (pdf) 


p0<s) 


exp 


When  the  signal  is  present  the  pdf  is 


pi(z) 


exp 


(2.40) 


(2.41) 
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In  (2.40)  on  represents  the  noise  amplitude  variance;  in  (2.41) 

2 u 

represents  the  signal-plus-noise  amplitude  variance. 

The  multidimensional  pdf  of  all  of  these  random  variables  (assumed 
to  be  independent)  is  given  by 

N-L  L 

P > *2  *****  * 71  Po^i>  71  P1(z.)  (2.42) 

i-1  j-1 

The  desired  probability  p is  given  by 
00 

P * f P(za)p0(z8)  dz8  (2.43) 

J0 

where  p(z  )is  the  probability  that  the  energy  sample  z in  a 
s s 

noise  slot  exceeds  the  energy  samples  in  all  the  other  slots. 

It  is  given  by 

p<v"  fo  1L(zs)  Fi<zs)  (2-44) 

where  FQ(zg)and  F^(zg)are  the  cumulative  distributions  correspond- 
ing to  the  pdf's  given  in  (2.40)  and  (2.41)  For  either  case  the 

(2.45) 


for  which  the  appropriate  subscript  can  be  added  to  o . 


distribution  has  the  generic  form 

F(z)  - / p(z)  dz 
*X) 


1 - exp 


2cr  ‘ 
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Using  these  results  in  (2*43)  we  have 


CO 

*0  l 


1 - exp 


2a; 


, N- 1-L 


1 - exp 


2a! 


exp 


2at 

0 

(2.46) 


De fining 


r “ an'Va 


2,  2 

0 /<ri 


(2.47) 


and  making  the  change  of  variables 


u - exp  { r 

( 2ctq 

we  obtain 


1 

P ■ f (1  - u)N"1"L(l  - ur)L  du 

*4 0 


Before  evaluating  this  integral  we  note  that 


(2.48) 


(2.49) 


crj2  » i GC2(0)  + Oq  (2.50) 

9 

In  this  expression  cr*  is  half  the  complex  noise  power,  C(0)  is 
the  peak  correlation  function  value,  and  G is  the  path  strength, 
i.e.,  G “ TiU2  where  is  the  complex  path  gain  of  the  itl1 
path.  Note  that  equal  strength  paths  have  been  assumed. 


dz 
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Recalling  that  the  noise  power  is  given  by 


- 2N0C<0) 

we  can  substitute  into  (2.47)  and  obtain 


(2.51) 


r - 


1 + 


L \TS 


LGCjOl  I 


2N, 


(2.52) 


where  is  the  duration  of  one  cycle  of  the  preamble.  It  is 
more  useful  to  write  (2.52)  in  the  form 


1 + 


Vo 


(2.53) 


where  y ^ the  power- to-noise  density  ratio. 

Returning  now  to  the  problem  of  evaluating  the  integral 
in  (2.49)  we  observe,  since  L cannot  be  a very  large  integer, 
that  there  is  no  computational  inconvenience  in  expanding  the 

v T 

term  (1-u  ) using  the  binomial  theorem.  We  obtain 


- £ O'*’1-1  / <•-■>"  •' 

j-0  Jo 


(L-j) 


du 


(2.54) 


where 


M - N - 1 - L 


The  mathematical  problem  becomes  one  of  evaluating  the  integral 
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(1  - x)  x dx 


This  is  accomplished  via  successive  integration  by  parts 
Initially  we  use  the  substitutions 


and 

dv  c xa  dx 
and  obtain 


I 

a 


(1-x) 


M x 


a+1 


a+1 


JO 


(-l)M  (1-x) 


M-l  x 


a+1 


a+1 


M 

a+1 


(1-x) 


M-l 


a+1 

x 


dx 


Proceeding  similarly  a second  time,  we  obtain 


(1  -x) 


M-2 


a+2 

x 


dx 


For  the  M 


th 


iteration, 


we  obtain  the  simple  result 


si 


Substituting  back  in  (2.54)  we  obtain 


M! 


p " E(j)  ST 

f|[r(L-j)  + il 


j-0 


1-1 


which  after  observing  that 


JJ  [r(L-j)+i]  - (M+l)!  JJ  [l+ 
i=l 


M+i  r 

77 

i-1 


and  using  the  fact  that 


M - N-l-L 


simplifies  to 


(2.61) 


(2.62) 


HI 


i. 


f ■ 
* 

l L 


I 


P = E (j)  M-L-r  1 -1 

j-o  77  1 + £iIrl1 


i-1 


(2.63) 


When  used  in  conjunction  with  (2.53)  and  (2.39),  this  result  yields 
the  rms  error  e^. 

For  a 31  length  sequence  and  a pulse  duration  of  1.39  ms, 
the  rms  frame  sync  error  is  presented  in  Figure  2.7.  The  error 
is  seen  to  decrease  rapidly  with  increasing  signal-to-noise- 
density  ratio  (SNDR)  and  number  of  paths.  When  the  SNDR  exceeds 
40  dB-Hz,  the  rms  error  remains  below  2 ms  in  all  cases. 
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Figure  2.7  RMS  Frame  Sync  Error  in  Excess  of 

Multipath  Error  for  an  L-Path  Channel 
with  Equal  Strength  Independently 
Fading  Paths 
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The  curves  presented  in  Figure  2.7  show  the  synchronization  j ; 

error  due  to  noise  only  and  should  be  thought  of  as  an  error  in  ] 

excess  of  the  multipath  error.  Inclusion  of  the  multipath  error  I 

] i 

requires  assumption  of  a specific  path  structure.  In  Figure  2.8  j ; 

i 1 

we  show  our  calculated  curves  for  a 2-path  channel  with  1 ms  \ \ 

between  the  paths  and  a 6-path  channel  with  0.5  ms  between  j 

adjacent  paths.  Plateaus  are  in  evidence  ahd  the  curves  are  seen  j 

to  be  more  tightly  clustered  together  than  in  the  previous  case.  i; 

j 

It  should  be  pointed  out  that  the  dependence  on  L in  j 

Figure  2.8  should  not  be  considered  as  indicative  of  the  dependence 

that  will  occur  in  practice.  The  wide  range  of  multipath  spreads  j \ 

: 1 

and  the  possibility  of  overlap  between  the  modes  allow  no  * j 

generalization  in  this  regard.  ) | 

j i 

■ i 

: ’ 

; j 

1 

: I 
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SECTION  3 


PREAMBLE  SIMULATIONS 


3.1  Signal  Presence  Detection  and  Doppler  Estimation 

3.1.1  Signal  Model 

Computer  sc.tware  for  simulating  modem  signal  presence  detection 
and  Doppler  estimation  has  been  developed.  A detailed  description  of 
the  modem  signal  processing  is  available  in  [3.1].  In  the  simulation 
these  operations  (bandpass  filtering,  frequency  shifting,  derivative 
estimation,  etc.)  are  carried  out  exactly  as  would  be  done  in  the 
actual  modem  although  floating  point  arithmetic  was  used.  In  the 
preamble  simulations,  the  channel  model  was  simplified  to  the  extent 
that  signal  structure  was  fixed  during  the  interval  (approximately 
350  ms)  required  for  these  operations . The  small  variation  which 
would  occur  in  a time-varying  channel  (bandwidth  on  the  order  of 
1 Hz)  can  be  expected  to  improve  rather  than  degrade  performance. 
Specifically,  if  s(t)  is  the  complex  envelope  which  would  be  re- 
ceived in  the  absence  of  distortion,  the  actual  received  signal  is 
of  the  form 


path 

r(t)  * 4^  aj  s(t_Tj)  + 


where  = delay  of  j 


th 


path 


(3.1) 


th 

relative  power  of  j path  satisfying 


(3.2) 
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During  signal  presence  detection/Doppler  estimation  the  signal 
consists  of  four  unmodulated  sinusoids  located  within  an  IF 
ban  dwidth  less  than  W(*2400)Hz,  as  shown  in  Figure  3.1: 


v— *\  j (2irf . t-Hp.  ) 

s(t)  - 72P  2^  e 


(3.3) 


Taken  together  with  (3.2)  this  defines  the  per  tone  received 
power  as  P.  Substituting  this  into  (3.1)  allows  us  to  write  the 
received  signal  in  the  simple  form 


<t}  ' X)  V 


j(2trfit4<pi) 


+ n(t) 


(3.4) 


where 


path 

:.«V2P  «y 

J-l 


(3.5) 


is  the  resultant  envelope  for  the  i tone.  In  an  exact  model, 
ou  and,  hence,  would  vary  slowly  during  the  signal  processing; 
this  effect  is  of  secondary  importance  and  has  been  omitted  in  the 
current  simulations. 

I.u  the  Vi..tulation,  r(t)  was  sampled  W times  per  second, 
corresponding  to  the  rate  at  which  complex  samples  are  processed 
in  the  actual  modem.  These  samples  are  thus  of  the  form 


Successive  noise  samples  n^j-j  arc  uncorrelated  zero-mean 
Gaussian  and  have  variance  given  by 


n(w)l 


2NqW 


(3.7) 


N 

where  __0  is  the  (two-sided)  power  spectral  density  of  the 

2 

noise.  Thus  the  SNR  of  the  waveform  being  processed  is  given  by 


N 


2P 


path 

E I-, 

J-A.. 


|n(w)l 


N0W 


(3.8) 


3.1.2  Variable  Bandwidth  Bandpass  Filters 

The  variable  bandwidth  filters  used  to  isolate  the  Doppler 
tones  consist  of  the  cascade  of  two  identical  first-order 
recursive  filters.  Consider  first  the  baseband  equivalent  of 
a single  filter: 

yk  * (1 ' 1*^-1  + ^ (3.9) 

The  transfer  function  of  this  filter  is  given  by 

- -ZI rr qai  (3-10> 

1 M - (H  - 1)  e JW1 

and 

, . l2  1 <3*n) 
|H^(jw)|  1 + y - y cos  o/r 

3-4 


where 


(3.12) 


-2  -1 
y m 2n  -2#i 

The  3 dB  frequency  of  the  cascade  may  be  found  by  setting 
|H.  (jte)  | 2 ■ ~Z  yielding 

1 J2 


w3dB 


or  equivalently 

7 * V2 - 1 
l-cos  u3dBr 

In  terms  of  y the  parameter  n is  given  by 


(3.13) 


(3.14) 


(3.15) 


The  resulting  amplitude  response  for  f3dB~75  Hz  is  shown  in 
Figure  3.2. 

In  the  simulation  an  equivalent  second-order  recursive 
filter  was  used  rather  then  separate  first-order  filters  and 
the  filtering  was  performed  at  bandpass.  For  an  assumed  center 
frequency  the  filtering  operation  is  defined  by 


Al6 


JV 


yk-l 


+ A2  e 


2jV 


yk-2  + A3*k 


(3.16) 


where 


A1  - 2(1-4)  - 2 


l _ /H:2v  -1 
y 


A2  - -(l-M)  - - 


1_yr+2v  -i 

y 


(3.17) 

(3.18) 


(3.19) 


To  achieve  a desired  3 dB  bandwidth  these  coefficients  must  be 
evaluated  using  Eq . (3.14). 


Filter  outputs  are  sampled  for  further  processing  at  an 

interval  sufficient  to  ensure  that  the  noise  in  successive 

samples  is  substantially  uncorrelated.  In  [3.1]  the  recommended 

separation  is  3t  where  t appears  in  the  equivalent  first-order 

c c 

continuous -time  filter  model  as 


H(f) 


j2.£+t^1 


(3.20) 


By  straightforward  manipulation,  it  can  be  shown  that  for  this 

model  t is  related  to  the  overall  3-dB  frequency  by 
c 


0.1024 
f3  dB 


(3.21) 
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j 3.1.3  Doppler  Estimation 

Although  the  Doppler  estimation  has  been  discussed  in 
* detail  elsewhere  [3.1],  we  present  a somewhat  different 

formulation  which  enables  a slight  reduction  in  the  number 
of  arithmetic  operations  required.  Consider  first  a single 
complex  envelope  sampled  at  time  kT; 

'*  zk  - \+5yk  (3.22) 

H * The  Doppler  estimate  based  on  this  signal  is 


A 


x,  y,  -y,  x, 
k^k  •'k  k 

2 2 

\ +yk 


(3.23) 


Si 


For  a constant-envelope  noise-free  signal  with  fixed  Doppler 
(z  = | z |eJa>r)  this  relationship  is  exact.  However,  the 
derivative  will  be  approximated  by  a simple  difference 

\ ’ T(veJ“T  2k-l)  (3’24) 


where  z,  1 represents  the  previous  sample  at  the  basic  sampling 
rate  and  the  factor  e-1  is  included  because  the  desired 
signal  imbedded  in  z(t)  is  the  output  of  a bandpass  filter 
centered  at  tu.  Substituting  (3.24)  into  (3.23)  we  find  that  the 
first  term  is  identically  zero,  yielding 


X Im  | e'JagV^_ij 
% " 2"T  |2k|2 
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(3.25) 


Calculation  of  the  numerator  of  (3 . 25) 'q  tires  6 multiplications 
and  3 additions,  while  the  original  approach  (calculating  x^ 
and  before  forming  requires  6 multiplications  and 

7 additions. 

At  N instants  in  time  (separated  by  approximately  t 
seconds),  the  numerator  and  denominator  of  (3.25)  are  computed 
for  each  of  the  four  bandpass  filter  outputs  and  accumulated 
separately.  Thus,  the  Doppler  estimate  may  be  regarded  as  a 
weighted  average  of  individual  Doppler  measurements 

N 4 


y y i z . i 2 v . 

Z— J L-U  1 n,  i 1 n,i 


n=l  i=l 
N 4 


(3.26) 


E E 'Zn,il 


n=l  i=l 


where  ^ represents  the  measurement  obtained  from  the  i 
tone  at  time  t=3ntc  and  ^ is  the  corresponding  complex 
envelope. 

3.1.4  Simulation  Results 

In  running  the  simulation  software,  the  following  parameter 
values  were  selected: 

Tone  Frequencies  (phases):  675.00(0),  1237.50(0), 

1800.0(0),  2362.50(180°) 

Noise  filter  center  frequencies:  956.25,  1518.75, 

2081.25,  2643.75 

Complex  samples  per  second:  3600 

Window  for  signal  presence  detection  ( t,  ] 

.0417  sec  (N-10)  V *IN' 

Signal  presence  threshold:  2.65 


it’l:  — ‘ - Jin'  L.,  !i 


Signal  presence  filter  bandwidth  = 75  Hz 


Doppler  estimation  filters: 

Stage  Bandwidth  (Hz) 

N 

Sample 

Separation  (ms) 

1 

75 

16 

4.17 

2 

16.2 

3 

18.89 

3 

3.5 

1 

87.78 

Duration  of  Doppler  estimation  (including 
transient  build-up):  0.3136  sec 

Initial  frequency  offset:  0 


The  parameters  most  important  to  the  design  are  the  bandwidths 
of  the  tone  filters  for  each  stage,  the  time  separation  between 
samples  taken  at  the  output  of  these  filters,  and  the  number  (N) 
of  samples  averaged.  Values  for  these  parameters,  as  shown 
above,  were  taken  from  the  preliminary  design  presented  in  [3.1]. 
The  sample  separation  must,  of  course,  be  an  integral  multiple  of 
the  basic  sampling  rate,  which  was  assumed  to  be  3600  Hz,  so  these 
values  differ  slightly  from  the  original. 

The  simulation  software  is  written  to  carry  out  signal 
presence  detection  and  Doppler  estimation  in  sequence.  However, 
if  signal  presence  detection  does  not  occur  during  the  first 
tWIN  seconc*s»  Doppler  estimation  is  not  attempted  and  the  trial 
is  recorded  as  an  acquisition  failure.  (In  the  actual  processor, 
acquisition  might  occur  subsequent  to  the  window,  but  Doppler 
estimation  would  then  overlap  the  wideband  sync  pulses  and  be 
substantially  degraded.)  Thus,  we  may  compare  the  observed 
probability  of  incorrect  dismissal  with  theoretical  calculations 
made  earlier.  (The  probability  of  false  alarm  was  also  predicted 
theoretically,  but  for  the  threshold  settings  of  interest  it  is 
too  small  to  observe  in  the  simulations.)  Recall  that  calculations 
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Figure  3.3  Simulated  Performance  of  Signal  Processing 

Detection  in  Flat  Fading  (Single-Path  Channel) 
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were  made  under  two  assumptions : completely  correlated  fading 
of  the  tones  and  independent  fading.  The  former  case  corresponds 
to  the  model  of  Eq.  0.1).  However,  as  the  number  of  paths  increases, 
the  tones  will  tend  to  become  uncorrelated „ Therefore,  we  have 
compared  simulation  results  for  the  two-  and  six-path  models 
with  the  predictions  based  on  independent  fading  in  Figure  3.4. 

These  results  also  show  good  agreement  in  the  region  for  which 
the  simulation  results  are  valid.  However,  the  effect  of 
correlation  between  the  tones  would  be  more  pronounced  at 
higher  values  of  SNDR. 

RMS  Doppler  estimation  error  as  a function  of  SNDR  is  shown 
in  Figures  3.5  and  3.6  for  the  two-  and  six-path  channel  models. 

Also  shown  are  theoretical  predictions  based  on  Eq.  (2.121)  of  [3.1]. 
In  general,  the  simulations  show  better  performance  than  the 
predictions.  This  is  because  Eq . (2. 121) assumes  that  in 
each  stage  the  tones  fall  at  the  3 dB  point  of  the  bandpass 
filter,  resulting  in  a loss  of  3 dB  in  SNR.  In  the  simulations, 
the  tones  were  located  at  the  center  of  the  first  stage  filters; 
in  subsequent  stages, tone  location  depends  on  the  estimate  made 
in  the  preceding  stage  during  that  particular  trial.  Thus,  the 
theoretical  predictions  are  overly  conservative  while  the 
simulations  more  accurately  represent  operational  conditions 
Examinations  of  the  curves  does  in  fact  show  a displacement 
of  roughly  3 dB  at  the  higher  SNDR's.  At  the  lower  SNDR's, 
performance  of  the  second  and  especially  the  third  stage  is 
degraded  as  a result  of  earlier  estimates  placing  the  signal, 
outside  of  the  filter  passband.  In  view  of  these  considerations, 
agreement  between  theory  and  simulation  seems  to  be  excellent. 
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PROBABILITY  OF  INCORRECT  DISMISSAL 


Figure  3.4  Simulated  Performance  of  Signal  Presence  Detection 
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3.2  Frame  Sync  Detection 


3.2.1  Wideband  Sync  Waveform 

The  wideband  sync  waveform  consists  of  a sequence  of 
bipolar  bandlimited  pulses  modulating  a carrier: 

v-1  j^Trf-t 

s(t)  »^2Pt  2 ymP(t‘^)e  (3.27) 

The  frequency  fg  is  chosen  so  that  the  pulse  spectrum  falls  within 
the  equipment  passband,  as  shown  in  Figure  3.7.  The  sequence  |y  j 
consists,  in  this  case,  of  three  repetitions  of  a length  31 
m-sequence,  chosen  for  its  desirable  autocorrelation  properties, 
as  discussed  in  Section.  2. 

In  the  receiver  signal,  processing  consists  of  an  approximate 
matched  filter  p(t),  followed  by  a tapped  delay  line  structure 
whlrh  forms  the  correlation  between  the  original  sequence 

and  appropriate  samples  of  the  output  of  the  matched  filter. 

This  structure  is  shown  in  Fig.  3.8  for  a pulse  separation  A 
equal  to  5 times  the  sampling  interval  t of  the  modem;  the 
choice  of  this  parameter  will  be  motivated  below.  The  output 
of  the  shift  register,  available  every  t seconds,  is  compared 
over  one  full  period  of  the  m-sequence  or  31  A seconds,  and  the 
largest  value  is  selected.  The  initiation  of  this  procedure  is 
timed  so  that  the  sync  waveform  is  present  at  the  input  during 
the  entire  interval  to  avoid  end  efforts.  This  is  the  basis  for 
the  receiver's  estimate  of  the  arrival  time  of  the  signal. 


0 


1 


Figure  3.7  Transmitted  Pulse  Spectrum  for 
Sync  Preamble  Pulse  Train 


Before  processing  the  sync  preamble  the  receiver  has  a 

timing  uncertainty  of  41.7  ms  resulting  from  the  time 

allowed  for  signal  presence  detection.  In  order  that  no 

timing  ambiguity  result  from  sync  preamble  detection,  it  is 

desirable  that  the  period  of  the  sync  preamble  be  somewhat 

greater  than  this  value.  For  a sampling  rate  of  3600  sec  \ 

the  31  bit  m-sequence  will  satisfy  this  goal  if  A > 4.84  t . 

s 

However,  A must  be  an  integral  multiple  of  t in  order  for  the 

s 

structure  of  Fig.  3.8  to  be  valid.  Thus,  we  are  led  to  the  value 

of  A ■ 5t  for  a preliminary  design.  Note,  however,  that  more 
s 

timing  accuracy  could  be  obtained  with  a shorter  pulse  width 

(A  ■ 4t  , for  example).  This  would  require  a longer  m-sequence 
s 

to  avoid  ambiguity  or  some  entirely  different  means  of  resolving 
the  resulting  ambiguity. 

The  transmitted  pulse  p(t)  was  selected  to  have  a band- 
limited  "raised-cosine"  spectrum  and  unit  power.  Thus  p(t) 

x ie 

and  its  transform  P(f)  have  the  form 


(3.28) 


(3.29) 


This  represents  a departure  from  [3.1]  where  it  was  proposed  that 
the  convolution  of  p(t)  and  the  receiver  pulse  matched  filter  have 
periodic  zeros  rather  than  p(t)  as  in  (3.28).  As  a consequence,  the 
sidelobe  levels  of  the  correlation  function  are  somewhat  degraded 
from  ideal.  However,  p(t)  and  the  matched  filter  are  more  readily 
realized  because  of  their  reduced  "tails".  Resources  available  did 
not  permit  a tradeoff  study  of  the  many  waveforms  available. 
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This  baseband  signal  is  shifted  so  that  it  occupies  a bandwidth 
of  2/ A located  within  the  passband  of  the  HF  SSB  equipment  as 

For  A 
t 


5t  as  established  above  and 
s 


shown  in  Fig.  3.7. 

36b 5 sec>  t*le  required  bandwidth  is  1440  Hz,  easily  within 


the  equipment's  capability.  In  selecting  the  center  frequency 
f^,  it  is  profitable  to  consider  the  receiver  signal  processing 
shown  in  Fig.  3.8.  Since  the  output  of  the  approximate  matched 
filter  will  be  centered  at  fg,  successive  samples  at  its  output 
will  be  shifted  in  phase  by  2fffQtg  radians,  and  successive 
samples  taken  from  the  shift  register  will  be  shifted  by  5*2frfQtg 
Combining  these  samples  in  the  summer  to  produce  a baseband 

output  would  ordinarily  require  multiplying  each  by  the  proper 

-j5  • 2*r  fnt 
J Os 


in  addition  to  multiplication  by  y 


m 


multiple  of  e 

However,  this. can  be  avoided  by  choosing  fjj  so  that  this  factor 
is  unity;  in  other  words 


5-27rf0ts  = k(2n) 


(3.30) 


where  k is  any  integer.  Selecting  k = 2 results  in  f^  * 1440  Hz 
so  that  the  transmitted  spectrum  falls  conveniently  between 
720  and  2160  Hz. 

Since  p(t)  in  Eq.  (3.28)  has  infinite  duration,  some  approxima- 
tion must  occur  in  the  receiver  pulse  matched  filter  in  Fig.  3.2. 
For  the  simulations,  this  filter  was  formed  by  truncating  p(t) 
to  9 samples  symmetrically  located  about  t * 0 and  using  these 
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samples  as  the  coefficients  of  a FIR  filter.  Since  the  next 
samples  in  either  direction  are  zero,  this  represents  an 
accurate  approximation  to  p(t).  Fig.  3.9  shews  the  output  of  the 
receiver  processing  (Fig.  3.8)  when  its  input  is  the  noise-free 
transmitted  signal. 

In  simulating  wideband  sync  detection,  a signal  model 
similar  to  that  used  in  simulating  signal  presence  and  Doppler 
estimation  was  used.  The  only  significant  difference  is  that 
the  transmitted  signal  structure  is  given  by  Eq.  (3.27)  rather 
than  by  Eq.  (3.3);  other  details  of  channel  structure  and  sampling 
are  identical. 

3.2.2  Simulation  Results 

Theoretical  predictions  of  sync  detection  performance  are  pre- 
sented in  Section  2.  To  compare  simulation  with  predictions,  we 
have  plotted  simulation  results  for  a one -path  channel  and  the  cor- 
responding (L « 1)  prediction  in  Fig.  3.10.  The  two  curves  are  in 
good  agreement  in  spite  of  the  different  assumptions  made  in  obtain- 
ing them. 

First,  the  theoretical  calculation  assumes  an  idealized  auto- 
correlation function  having  zero  width,  while  the  simulation  uses 
the  actual  autocorrelation  function  for  the  bandlimited  transmitted 
waveform,  including  the  effect  of  the  approximate  matched  filter  in 
the  receiver. 

Second,  simulation  of  sync  detection  is  performed  in  sequence 
following  signal  presence  detection  and  Doppler  estimation.  Any 
error  in  the  final  Doppler  estimate  is  present  during  sync  detection, 
an  effect  not  considered  in  the  theoretical  predictions. 
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Figure  3.9  Discrete  "Autocorrelation"  Resulting  from 
Approximate  Matched  Filter 


A third  effect  occurs  when  more  than  one  path  is  present. 

With  multipath  spreads  of  less  than  3 ms,  interference  (construc- 
tive or  destructive)  between  paths  will  occur  in  the  receiver 
processing.  This  effect  is,  of  course,  modeled  in  the  simulation 
but  not  taken  into  account  in  the  theoretical  predictions.  For 
this  reason,  we  will  not  compare  theory  and  prediction  when  mul- 
tiple paths  are  present. 

Simulation  results  for  the  two-path  and  six-path  channel  models 
are  plotted  on  Fig.  3.11.  When  additional  paths  are  included,  three 
effects  occur.  First,  the  multipath  spread  of  the  channel  itself 
causes  a residual  error  in  the  sync  estimation  which  will  be  present 
even  at  high  SNR.  This  effect  is  evident  in  the  two-  and  six-path 
curves  of  Fig.  3.11.  Second,  the  presence  of  additional  independent 
paths  adds  diversity  which  tends  to  reduce  the  probability  of  detec- 
tion errors  due  to  noise  peaks.  However,  this  effect  is  offset  to 
some  degree  by  the  effective  loss  of  power  resulting  from  the  fact 
that  the  receiver  partially  resolves  the  paths.  It  turns  out  that, 
for  rms  errors  larger  than  the  residual,  performance  in  the  two- 
path  case  is  roughly  3 dB  worse  than  in  the  one -path  situation. 

For  the  six-path  channel,  the  paths  are  spread  uniformly  over  a 
2.5-ms  delay  interval.  The  width  of  the  autocorrelation  function 
is  such  that  six  individual  paths  cannot  be  resolved  by  the  receiver. 
However,  received  energy  is  smeared  over  this  interval,  resulting  in 
further  degradation  from  the  two-path  case,  as  shown  in  Fig.  3.11. 
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3.3  Doppler  Tracking 

3.3,1  Decision-Directed  Doppler  Tracking 

Computer  simulations  of  Doppler  tracking  in  the  AN DVT  modem 
were  performed.  The  goals  of  this  effort  were  to  evaluate  the 
suitability  of  the  decision-directed  error  signal  proposed  previously 
[3.1]  and  to  design  a suitable  loop  filter.  A block  diagram 
of  the  tracking  loop  is  shown  in  Fig.  3.12a.  At  time  pT, the  Doppler 
srror  signal  € is  derived  by  the  demodulator  as  a sum  of  terms 
obtained  from  each  of  the  39  demodulated  tones: 


€ 


P 


k=l 


(3.31) 


If  S denotes  the  complex  output  of  the  k'th  tone  filter  at 
kp 

time  pT,  the  DPSK  demodulation  process  involves  forming  the 
product 


S*  S , = R.  + iQ,  tJ.oz 

bkp  k,p-l  Kp  Jykp 

The  signs  of  the  real  and  imaginary  parts  and  Q^p  give  the 
decisions  on  the  transmitted  data.  In  terms  of  these  quantities, 
the  Doppler  error  signal  for  this  tone  is  given  by 

fkP  ■ <'V  ■ % u sgn  (V sga  (V  (3’33 

It  has  been  shown  earlier  [3.1]  that  with  slow  fading  and  in  the 
absence  of  noise 


■IrliJtfcvAV.. u.: . 
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cfc  = Jl  |T(kF,pT)|  iT[kF,(p-l)T]|  sin  2rr(i'0  - i^)  (3.34) 

where  T(f,t)  is  the  time-varying  transfer  function  of  the  channel. 

The  design  goal  has  been  operation  of  the  tracking  loop  in 
the  presence  of  a 3.5  Hz/sec  frequency  offset.  A second-order  loop 
filter  (see  Fig.  3.12b)  was  selected  to  give  zero  steady-state 
tracking  error  for  a ramp  input.  Additional  goals  which  have  been 
satisfied  are: 

(1)  Successful  tracking  at  the  initiation  of  a 

- 3.5  Hz/sec  ramp  in  frequency 

(2)  Successful  tracking  at  the  initiation  of  a 
1 Hz/sec  step  in  frequency 

Figure  3.13  shows  acquisition  of  a 3.5  Hz/sec  ramp  (for  a 
noise-free,  nonfading  signal)  for  loop  constants  K=0.05  and 
a*  0.03.  Figure  3.14  shows  the  acquisition  of  a 1-Hz  frequency 
step  under  these  same  conditions.  These  values  of  loop  constants 
a and  K represent  a compromise;  increasing  loop  bandwidth  would 
improve  transient  response  but  increase  susceptibility  to  noise. 

.Since  Doppler  tracking  is  a sampled  data  process  operating 
at  the  frame  rate,  the  channel  model  used  for  simulating  signal 
presence  detection,  Doppler  estimation,  and  sync  detection  would 
be  inappropriate.  The  channel  model  used  for  simulating  long-term 
Doppler  tracking  has  been  described  earlier  [3.1,  Appendix  B] . 

The  only  modification  necessary  has  been  the  inclusion  of  the 
frequency  ramp  representing  the  Doppler  offset. 
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FRAMES 


Figure  3.13  Response  to  a Ramp  in  Frequency 
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3.3.2  Simulation  Results 


I 

I 

To  provide  a comparison  for  the  decision- feedback  tracking 
J simulation  results,  some  simulations  of  conventional  single-tone 

tracking  have  been  carried  out.  In  this  scheme  an  unmodulated 
tone  of  equal  amplitude  is  added  to  the  ANDVT  format  at  562.5  Hz. 

At  the  receiver  this  tone  is  selected  by  a bandpass  filter  and 
passed  through  a discriminator  to  generate  the  error  signal  for 
the  Doppler  tracking  loop.  The  structure  and  bandwidth  of  this 
loop  were  the  same  as  that  used  for  decision  feedback  tracking 
so  that  the  single-tone  tracking  simulation  differs  only  in  the 
error  signal  used.  In  the  single-tone  tracking  scheme  the  bandwidth 
of  the  bandpass  filter  together  with  the  received  noise  level 
establish  the  noise  power  at  the  input  to  the  discriminator.  As 
signal  fading  occurs,  the  instantaneous  signal-to-noise  ratio  of 
the  discriminator  will  fall  below  threshold,  causing  an  error 
signal  consisting  primarily  of  noise  to  be  supplied  to  the  tracking 
loop  until  the  received  signal  level  increases  again.  Clearly, 
the  bandpass  filter  should  be  made  as  narrow  as  possible  to 
L minimize  this  effect. 

jj  We  first  consider  simulation  results  for  the  two-path  chan- 

nel. Shown  in  Fig.  3.15  are  rms  tracking  errors  as  a function 
I'  of  E,/Nn  for  bandwidths  of  + 22.22  and  + 44.44  Hz  and  for  the 

decision-directed  scheme.  Table  3-1  shows  the  values  of  E,  /Nrt 

b o 

below  which  loss  of  lock  occurs  with  each  technique.  Together 
with  Fig.  3.15,  these  results  indicate  that,  of  the  three  systems 
being  compared,  decision-directed  tracking  gives  the  smallest 
tracking  error,  while  single- tone  tracking  with  a bandwidth  of 
+ 22.22  Hz  gives  the  lowest  threshold. 
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Figure  3.15  Tracking  Performance  of  Various  Doppler  Correction  Loops 
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The  effect  of  these  tracking  loops  on  error  rate  is  shown 
in  Fig.  3.16.  Performance  of  the  decision-directed  technique  falls 
on  top  of  the  reference  curve  at  all  E^/Nq  for  which  the  loop  is 
in  lock.  The  same  is  true  of  single-tone/discriminator  tracking 
with  a bandwidth  of  + 22.22  Hz.  For  a bandwidth  of  + 44.44  Hz 
the  curve  shown  in  Fig.  3.16  was  obtained,  showing  a performance 
loss  of  1 - 2 dB.  It  should  be  noted  that  in  the  latter  simulation 
for  the  single-tone/discriminator  no  correlation  was  assumed 
between  fading  of  the  Doppler  tone  and  fading  of  the  data  tones. 

The  effect  of  this  assumption  on  the  results  is  expected  to  be 
small. 

We  next  consider  simulation  results  for  the  six-path  channel. 
Shown  in  Fig.  3.17  are  rms  tracking  errors  as  a function  of  Eb/NQ 
for  bandwidths  of  + 22.22  and  + 44.44  Hz  and  for  the  decision- 
directed  scheme.  Table  3-2  shows  the  limiting  values  of  E^/Nq 
below  which  successful  tracking  is  not  possible.  Together  with 
Fig.  3.17  these  results  indicate  that,  for  the  six-path  channel, 
decision-directed  tracking  is  superior  to  single-tone/discriminator 
tracking  in  both  threshold  and  rms  tracking  error. 

The  effect  of  these  tracking  loops  on  error  rate  is  shown 
in  Fig.  3.18.  The  curve  for  discriminator  tracking  with  a bandwidth 
of  + 22.22  Hz  falls  on  top  of  the  reference  curve  at  all  Eb/Ng 
for  which  the  loop  is  in  lock.  The  same  is  true  of  decision- 
directed  tracking  except  for  values  of  E,  /Nft  just  above  threshold. 
With  discriminator  tracking  in  a + 44.44  Hz  bandwidth,  there  is 
significant  degradation  in  error  rate  at  all  E^/Ng. 
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Bandwidth  ■ 1 Hz 

Loop  Constants:  K ■ 0.05,  a ■ 0.03 

Figure  3.16  Simulated  Raw  Bit  Error  Rates  showing  Effect  of 
Closed  Loop  Frequency  Tracking 
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Figure  3.17  Tracking  Performance  of  Various  Doppler  Correction  Loops 
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Figure  3.18  Simulated  Raw  Bit  Error  Rates  Showing  Effect  of 
Closed  Loop  Frequency  Tracking 


[3.1]  D.  Chase,  et  al.,  "HF  Modem  Design  for  the  Advanced 
Narrowband  Digital  Voice  Terminal  (ANDVT),"  November 
1977,  Final  Report  on  Contract  N00039-77-C-0190. 


SECTION  4 

CODING  AND  MODULATION  FORMAT  FOR  THE 
KG  SYNCHRONIZATION  SEQUENCE 

4.1  Introduction 

The  results  obtained  during  the  previous  effort  indicated  that 
the  (250,100;46),  (252, 128;37) , and  (250,150;27)  BCH  codes,  combined 
with  eighth-order  diversity,  represent  attractive  candidates  for 
providing  the  COMSEC  synchronization.  In  addition  to  these  long 
BCH  codes,  it  is  of  interest  to  consider  the  use  of  several  (24, 12; 8) 
Golay  codewords  to  protect  the  KG  synchronization  sequence.  For 
example,  if  ten  Golay  codewords  are  used  with  eighth-order  diversity, 
one  may  view  the  overall  code  as  a (240,120),  which  may  be  compared 
to  long  BCH  codes.  In  this  section  we  first  compare  these  coding 
approaches  by  using  analysis  techniques  which  are  based  on  a Rayleigh 
fading  channel  model  with  code  symbols  assumed  to  be  independent. 
These  theoretical  results  are  shown  to  be  somewhat  optimistic  when 
compared  to  simulated  results  obtained  for  the  two-path  and  six- 
path  HF  models  used  throughout  this  project.  In  fact,  the  theore- 
tical results,  given  in  Section  4.2,  indicate  that  a long  BCH  code- 
word with  just  hard  decoding  is  somewhat  superior  to  the  use  of 
multiple  Golay  codes  with  soft  decoding.  The  simulation  results, 
given  in  Section  4.3,  which  include  the  correlation  between  code  sym- 
bols, renders  the  multiple  Golay  codes  with  soft  decoding  superior 
to  the  long  BCH  codes  decoded  with  just  binary  information.  The 
recommended  coding  techniques  for  th^  KG  synchronization  sequence 


is  also  given  in  Section  4.3.  It  should  be  pointed  out  at  the 
onset  of  this  section  that  both  the  multiple  Golay  coding  approach 


i 


and  the  long  BCH  codes  are  attractive  choices  since  these  codes 
provide  performance  capabilities  such  that  KG  synchronization 
can  be  achieved  even  when  the  channel  conditions  are  such  that 
the  quality  of  the  digitized  voice  would  be  unacceptable. 

4.2  Theoretical  Evaluations  Based  on  the  Independent 

Rayleigh  Fading  Model 

We  have  begun  this  investigation  by  obtaining  the  performance 
of  four  candidate  KG  codes  for  the  COMSEC  synchronization.  Any 
remaining  code  bits  with  a code  block  may  be  used  for  net  control. 

The  codes  are  a BCH  (250,100;46)  code,  a BCH  (252,128;37)  code,  a 
BCH  (250, 150; 27)  code,  and  ten  (24, 12; 8)  Golay  codewords,  which  can 
be  viewed  as  a (240,120)  code. 

Both  hard  decoding  and  soft  decoding  are  investigated.  The 
results  are  compared  on  the  basis  of  both  decoded  bit  error  rate  and 
block  error  rate.  In  addition,  both  slow  and  fast  fading  cases  are 
considered.  The  results  are  plotted  in  Figures  4.1  through  4.8. 

In  addition  to  comparing  performance,  it  is  of  interest  to 
modify  the  processing  time  study  for  Contract  N00039-77-C-0190  to 
provide  results  for  the  (252,128)  code  and  (250,150)  code  as  well 
as  results  already  obtained  for  the  (250,100)  code.  The  time  re- 
quired for  hard-decision  decoding  of  the  (24,12)  Golay  code  has  also 
been  isolated  from  the  study  of  the  processing  time  for  a 16-iteration 
soft-decision  decoding  of  the  (24,12)  extended  Golay  code. 

Table  4-1  summarizes  the  processing  times  for  the  various  codas 
and  algorithms.  The  processing  time  is  given  in  terms  of  the  t^y 

and  t parameters  which  have  been  described  in  the  previous 

memory 

final  report.  A value  of  350  ns  is  assumed  in  Table  4-1  to  provide 
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TGTPL  SNR  PER  RECEIVED  TONE  CET/NOl  IN  DB 

Figure  4.1  Hard-Decoded  Bit  Error  Rates  of  KG  Codes  over  a Slowly  Fading 
HF  Channel  with  4-Phase  DPSK  Modulation  using  Eighth-Order 
Post-Detection  Equal  Gain  Diversity  Combining 


Figure  4.2  Hard-Decoded  Block  Error  Rates  of  KG  Codes  over  a Slowly 


Fading  HF  Channel  with  4-Phase  DPSK  Modulation  using 
Eighth-Order  Post-Detection  Equal  Gain  Diversity  Combining 
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Figure  4.3 


TQTPL  SNR  PER  RECEIVED  TONE  CET/N01  IN  OB 

Estimated  Soft-Decoded  Bit  Error  Rates  of  KG  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation  using 
Eighth-Order  Post-Detection  Equal  Gain  Diversity  Combining 
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TOTAL  SNR  PER  RECEIVED  TONE  CET/NO)  IN  DB 


Figure  4. 


4 Estimated  Soft-Decoded  Block  Error  Rates  of  KG  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation  using 
Eighth-Order  Post-Detection  Equal  Gain  Diversity  Combining 
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TOTRL  SNR  PER  RECEIVED  TONE  [ET/NG)  IN  DB 

Figure  4.5  Hard-Decoded  Bit  Error  Rates  of  KG  Codes  over  a Time-Varying 
HF  Channel  with  4-Phase  DPSK  Modulation  (Frame  Rate  “ 44.44 
Frames /Second,  B “ 5 Hz)  using  Eighth-Order  Post-Detection 
Equal  Gain  Diversity  Combining 
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TOTPL  SNR  PER  RECEIVED  TONE  CET/NQ]  IN  DB 

Figure  4.7  Estimated  Soft-Decoded  Bit  Error  Rates  of  KG  Codes  over  a 

Time-Varying  HF  Channel  ^ith  4-Phase  DPSK  Modulation  (Frame 
Rate  “ 44.44  Frames /Second , B » 5 Hz)  using  Eighth-Order 
Post-Detection  Equal  Gain  Diversity  Combining 
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Block  Error  Rate 


TOTAL  SNR  PER  RECEIVED  TONE  CET /NO]  IN  DB 


Figure  A. 8 Estimated  Soft-Decoded  Block  Error  Rates  of  KG  Codes  over  a 
Time -Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  - 44.44  Frames /Second,  B-5  Hz)  using  Eighth-Order 
Post-Detection  Equal  Gain  Diversity  Combining 
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TABLE  4-1 


DECODER  PROCESSING  TIMES 


Code  and  Algorithm 

Processing  Time 

Estimate  for 

£ * ‘alu  * 350  08 

(249, 101; 45)  BCH 
(Hard  Decode) 

104837  tAT„  + 36136  t 

ALU  memory 

49.34  ms 

(252,128;37)  BCH 
(Hard  Decode) 

84292  tATTt  + 28678  t 

ALU  memory 

39.54  ms 

(250,150;27)  BCH 
(Hard  Decode) 

58504  tATtt  + 19546  t 

ALU  memory 

27.32  ms 

10  Codewords  using 
a (24,12)  Extended 
Golay  Code 
(16  Iteration  Soft 
Decode) 

134560  tATJT  + 28500  t 

ALU  memory 

56.13  ms 

10  Codewords  using 
a (24,12)  Extended 
Golay  Code 
(Hard  Decode) 

6480  tATIT  + 530  t 

ALU  memory 

2.45  ms 

an  approximate  timing  estimate.  These  processing  times  also  favor 
the  use  of  a long  BCH  code  with  hard  decoding,  insread  of  the  use 
of  a multiple  soft-decoded  Golf.y  coding  approach. 


However,  before  a firm  recommendation  can  be  made,  a complete 
simulation  over  realistic  channel  models  is  necessary.  These 
results  are  given  in  the  following  section. 

4.3  Simulation  Results  for  the  Two-Path  and  Six-Path 
Channel  Model 


In  this  section  a detailed  simulation  is  conducted  for  the 
(252,128;37)  BCH  code  and  11  (24,12;8)  Golay  codes.  Thus,  the  BCH 
code  protects  128  bits  which  may  be  used  of  KG  synchronization  net 
control,  and  the  11  Golay  codes  protect  132  bits. 


Transmission  of  the  encoded  sequence  is  to  be  accomplished  by 
a 39  parallel  tone  modem  using  4-phase  DPSK  modulation.  Seventy- 
eight  bits  of  data  are  relayed  during  each  22. 5 -ms  frame.  Since 
the  block  length  of  the  encoded  sequence,  N,  must  be  greater  or 
equal  to  128,  more  than  one  frame  will  be  required  to  transmit  each 
data  sequence.  For  the  purposes  of  providing  time  diversity,  each 
encoded  sequence  will  be  repeated  L times.  The  number  or  frames 
required  to  transmit  the  multiple  repeats  of  the  encoded  sequence 
is  simply 
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where  f,"j  denotes  the  ceiling  operator.  The  re- 
ceiver uses  equal  gain  predetection  combining  to  assemble  the 
results  of  the  L sequence  transmissions. 


A block  code  selected  for  the  encoding  process  must  satisfy 
the  dual  criteria  of  excellent  performance  at  the  signal-to-noise 
ratios  of  interest  and  relatively  simple  encoding  and  decoding 
algorithms . 


A simulation  of  this  type  slightly  favors  the  BCH  code  since 
fewer  bits  are  protected.  However,  it  is  expected  that  the 
results  for  the  Golay  codes,  when  only  120  bits  are  protected, 
will  be  quite  similar  to  these  simulated  results. 

4-J.2 


One  of  the  protection  schemes  of  interest  incorporates  a 
(252,128;37)  BCH  code.  The  bits  representing  the  first  trans- 
mission of  the  encoded  KG  synchronization  sequences  are  placed 
on  the  tone  array  as  diagrammed  in  Figure  4.9.  All  subsequent 
repetitions  of  the  KG  codeword  are  placed  on  the  tone  array  in 
packed  format.  An  example  of  the  packing  for  the  second  KG  code- 
word occurs  during  frame  i+3  in  Figure  4.9. 

f.L 

Each  of  the  L repetitions  of  the  j (lsjsN)  codeword  bit 
will  be  combined  before  demodulation.  In  order  to  benefit  from  the 
available  inband  diversity,  the  L bit  repetitions  should  appear  to 
be  evenly  dispersed  across  the  tone  array.  An  example  of  this 
dispersal  for  the  mcdulation  format  of  Figure  4.9  is  depicted  in 
Figure  4.10.  Tone  placements  of  the  first  codeword  bit  for  a maxi- 
mum of  eight  repetitions  are  easily  discerned.  Tone  positions  of 

i.L 

the  j bit  can  be  easily  determined  by  right-circularly-shifting 
the  tone  array  index  in  Figure  4.10  by  j-1. 

An  alternate  protection  scheme  requires  the  use  of  11  (24, 12; 8) 
Golay  codewords  to  form  a single  repetition  of  an  encoded  KG  syn- 
chronization sequence  data  block.  The  maximum  number  of  data  bits 
which  can  be  protected  is  132;  the  total  duration  of  the  encoded 
sequence  is  264  bits.  The  modulation  format  bit  assignments  for 
this  scheme  are  presented  in  Figure  4.11. 

The  11  Golay  codewords  comprising  each  KG  codeword  are  inter- 
leaved. However,  repetitions  of  each  KG  codeword  are  packed  into 
the  done  array  in  a manner  similar  to  that  described  for  the  BCH 
protection  technique.  This  results  in  the  j bid  of  each  of  the 
L KG  codeword  repetitions  being  placed  on  tones  which  are  separated 
in  frequency.  In  Figure  4.12,  tone  placements  of  the  L repetitions 
(LS8)  of  the  first  bit  of  the  first  Golay  codeword  is  indicated. 
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Modulation  Format  Bit  Assignments  for  KG  Synchronization  Sequence  when 
(252, 128; 37)  BCH  Code  Is  Employed 


Figure  4.12  The  Index  of  the  Tone  on  which  the  First  Bit  of  8 Consecutive 
Synchronization  Sequence  Codewords  llx(24,12;8)  Is  Placed 


The  tone  placement  of  the  24  bits  of  the  first  Golay  codeword  are 
depicted  in  Figure  4.13.  Both  diagrams  demonstrate  that  acceptable 
tone  dispersement  of  related  bits  has  been  achieved. 

To  demonstrate  the  relative  robustness  of  the  two  encoding 
techniques,  software  has  been  written  to  simulate  the  transmission 
of  the  KG  synchronization  sequence  over  an  HF  channel.  Two  chan- 
nel models  will  be  used  for  evaluation  purposes: 

Model  1 - The  HF  channel  is  represented  by  two  equi- 

strength  paths  separated  by  1 ms.  Each  path 
is  assumed  to  have  a double-sided  Doppler  band- 
width of  1 Hz.  This  model  is  representative 
of  channel  behavior  when  the  modem  is  operating 
near  the  maximum  usable  frequency  (MUF) . 

Model  2 - The  HF  channel  is  represented  by  a 6 -path 
equi-strength  split  ray  model  with  a total 
multipath  spread  of  2.6  ms.  Each  path  is 
assumed  to  be  resolvable  into  two  components 
separated  by  0.1  ms.  The  double -sided  Doppler 
bandwidth  per  path  is  1.5  Hz.  This  model  is 
representative  of  channel  behavior  when  the 
modem  is  operating  significantly  below  the 
maximum  usable  frequency;. 

A summary  of  the  details  of  both  channel  models  appears  in  Table  4-2. 

The  result  of  the  simulations  is  a statistical  summary  of  the 
decoded  KG  synchronization  sequence  bit  and  block  error  rates.  Of 
course,  the  block  error  rates  are  computed  on  the  basis  of  the  entire 
KG  codeword.  For  the  case  of  the  Golay  protection  scheme,  a KG  block 
error  occurs  whenever  one  or  more ‘'of  the  11  Golay  codewords  is  de- 
coded incorrectly.  For  the  case  of  the  BCH  protection  scheme,  a KG 
codeword  is  identical  to  the  block  error  rate  for  the  BCH  codeword. 
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TABLE  4-2 

SUMMARY  OF  THE  DELAY  AND  DOPPLER  CHARACTERISTICS 
FOR  TWO  SIMULATED  HF  CHANNEL  MODELS 


Statistics  have  bean  kept  for  both  hard  and  soft  decoding  of 
the  Golay  codewords.  The  soft  decoding  algorithm  employed  is 
Chase's  Algorithm  2.  This  method  requires  a maximum  of  16  hard 
decoder  Iterations  to  complete  a soft  decision  decode  of  the 
(24, 12; 8)  Golay  codeword.  Unfortunately,  the  large  minimum  distance 
of  the  (252, 128; 3 7)  BCH  code  renders  impractical  any  contemplated 
application  of  Chase's  Algorithm  2.  As  a result,  only  a hard 
decoding  simulation  has  been  done  for  the  BCH  protection  scheme. 

In  Figures  4.14  through  4.19,  block  error  rates  as  a function 
of  signal-to-noise  ratio  per  received  tone  (E^/Nq)  are  presented. 

The  curves  are  drawn  with  parametric  dependence  on  encoding  tech- 
nique, HF  channel  model,  and  the  number  of  KG  codeword  repetitions. 
Figures  4.20  through  4.25  present  bit  error  rates  as  a function  of 
signal-to-noise  ratio  per  received  tone  (Et/Nq)  with  similar 
parameterizations . 

Both  the  BCH  technique  and  the  Golay  technique  demonstrate 
approximately  equal  performance  when  the  Golay  codewords  are  hard 
decoded.  This  is  especially  true  for  the  two-path  HF  channel 
model.  However,  the  BCH  technique  does  suffer  an  increased  disad- 
vantage with  respect  to  the  hard-decoded  Golay  technique  for  the 
six-path  split  ray  HF  channel  model. 

As  noted  earlier,  the  relationship  between  the  three  SNR's 
used  in  this  report  is  given  by 


4? 

A simple  two  iteration  soft  decoding  algorithm,  based  on  one  con- 
ventional hard  decode  and  one  additional  hard  decode  with  the 
37  leas;  reliable  bits  complemented,  may  prove  useful  if  this 
code  is  finally  selected  for  the  ANDVT. 
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Figure  4.14  Block  Error  Rate  for  the  (252, 128; 37)  BCH  Code  Assuming  a 
2-Path  HF  Channel  Model  (L  ^ No.  of  KG  Codeword 
Repetitions)  a. 22 


N 


( 

Figure  4.15  Block  Error  Rate  for  a Hard-Decoded  Golay  Protection 
Scheme  Assuming  a 2-Path  HF  Channel  Model  (L  No. 
of  KG  Codeword  Repetitions) 

4-23 


• I * I I 

»-3  M ! ‘ > 


nr  j ! i ! Mil  i nr  j = s \ 


: > I l : I i ; l i : 


MiMji!  u 


. * . i 


* i ' • > * 


! : i • 1-fl 

i i i III 

-h-r--nr 


r H ! : ! 


1 3 x 10' 


LLii! J. 

3 -2 


! ! j |j  | i ; j f 

t]  I j ! i 

-i-i.  . : i.:  i 1 l : . 


nr-\| 


ET/N0 


ilium:  i. 

4 5 


6 ' 


Figure  4.16  Block  Error  Rate  for  a Soft-Decoded  Golay  Protection 
Scheme  Assuming  a 2-Path  HF  Channel  Model  (L  No. 
of  KG  Codeword  Repetitions)  ^-24 
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Figure  4.18  Block  Error  Rate  for  a Hard-Decoded  Golay  Protection 
Scheme  Assuming  a 6 -Path  Split  Ray  HF  Channel  Model 
(L  ■>  No.  of  KG  C^c^grord  Repetitions) 
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Figure  4.19 


Block  Error  Rate  for  a Soft-Decoded  Golay  Protection 
Scheme  Assuming  a 6-Path  Split  Ray  HF  Channel  Model 
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Codew^r^^Repetitions) 


Figure  4.21  Bit  Error  Rate  for  a Hard-Decoded  Golay  Protection 
Scheme  Assuming  a 2-Path  HF  Channel  Model  (L  ■>  No 
of  KG  Codeword  Repetitions) 


Figure  4.23  Bit  Error  Rate  for  the  (252,128;37)  BCH  Code  Assuming 
a 6-Path  Split  Ray  HF  Channel  Model  (L  ■>  No.  of  KG 
Codeword  Repetitions)  4-31 
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Figure  4.24  Bit  Error  Rate  for  a Hard -Decoded  Golay  Protection 
Scheme  Assuming  a 6-Path  Split  Ray  HF  Channel  Model 
(L  No.  of  KG  Cgd^yord  Repetitions) 


+ 33.8 


and,  thus, 


For  acceptable  quality  speech  on  a two-path  model,  anE^/Ng  of  14  dB 
is  required,  and  an  E^/N^  of  11  dB  is  required  for  the  six-path 
model.  Thus,  both  the  long  BCH  code  and  the  soft-decoded  Golay 
code,  used  with  eighth-order  diversity,  represent  attractive  codes 
for  KG  synchronization.  However,  the  Golay  code  with  soft  decoding 
offers  somewhat  better  performance  than  that  obtained  for  the  long 
BCH  code  with  just  hard  decoding.  Based  on  performance  results 
alone,  the  Golay  coding  approach  is  recommended.  A choice  of  ten 
Golay  codes,  protecting  120  bits,  appears  somewhat  more  suitable 
than  the  eleven  Golay  codes  actually  simulated.  The  computation 
is  somewhat  greater  for  the  Golay  coding  approach,  as  indicated  in 
Table  4-1 J however,  this  increase  in  computation  must  be  weighted 
against  the  additional  software  program  required  for  decoding  the 
long  BCH  code. 

It  shduld  be  noted  that  the  early  analysis,  assuming  independent 
Rayleigh  fading,  indicated  that  the  BCH  technique  would  be  superior 
to  the  soft-decoded  Golay  scheme.  The  importance  of  detailed  simu- 
lations such  as  undertaken  above  should  be  evident. 
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APPENDIX  A 


APERIODIC  CORRELATION  FUNCTION 


In  this  appendix  we  briefly  summarize  the  development  of 
the  aperiodic  autocorrelation  function  and  record  for  future 
reference  (and  comparison  with  more  general  results)  the  relevant 
expressions. 

The  transmitted  preamble  is  the  result  of  convolving  an 
impulse  train  of  length  K-l  with  a pulse  shaping  filter  p(t)  to 
obtain 

K-l 

z(t)  - E y p(t-kT)  (1) 

k-0  K 

At  the  receiver  the  preamble  is  processed  with  a pulse-shaping 
filter  q(-t)  and  a discrete  matched  filter.  The  output  of  the 
receiver  pulse-shaping  filter  is  given  by 

K-l 

Vt)  - E y.r(t-kT)  (2) 

1 k-0  K 


where 


r(t>  " J q(T)p(t  + T)  dT  (3) 

The  output  w(t)  of  the  discrete  matched  filter  is  the 
result  of  convolving  the  impulse  response  of  the  discrete  matched 
filter  with  (2) . We  have 


A-l 


*ytr*»mv.r 


W(t) 


'll  t 

^i-(K-l)  -1  llj-0  J 


r (t-t-jT)  dr 


" a 

£ E y^y^ft  - <i+j)T] 

i—  (K-l)  j-0  1 J 


K-l  K-l 
E E 
4-0  j-0 


• ^ ^ 

E y*y.r[t  - (j-4)Tl 
j-0  * J 


To  simplify  we  use  the  function  definition 


[a,b] (m) 


Then  w(t)  becomes 


+®  +«  | 

w(t)  - £ £ y*y  r[t 

1 — — 00  j =t»OD  J 


a < m < b 


otherwise 


- (j-i)TJ[0,  K-l] (1) [0,  K-l] (j) 


which,  with  a change  of  variables  k-j-1,  becomes 

+®°  . 

w(t)  - £ C(k)r(t-kT) 


where 


C(k)  - E Ytyi+k[0.  K-U(i)tO,  K-ll  <1  +k) 


E yty14k[0,  K-U(l)t-k,  K-l-k|(l)  (8) 


from  which  we  finally  obtain 


c(k) 


K-l-k 

A 


K-l 

{ t *k 


0 < k < K-l 


-(K-l)  £ k < 0 (9) 


Ik  | > K-l 


The  above  expression  can  be  simply  manipulated  to  obtain 
the  general  relation 


C(k)  - C(-k) 


(10) 


which  can  be  used  to  write  (9)  in  alternative  forms. 

Note  that  (9)  is  the  output  of  a non-realizable  filter. 
This  is  because  the  matched  filter  used  in  (4)  was  noncausal. 
To  analyze  practical  systems  so  that  we  can  keep  track  of  the 
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timing  it  is  necessary  co  insert  a delay  in  the  matched  filter 
impulse  response.  A shift  of  (K-l)T  puts  the  first  non-zero 
value  of  the  matched  filter  impulse  response  at  k*0.  C(k)  is 
affected  similarly.  Instead  of  running  from  -(K-l)  to  +(K-1) 
the  matched  filter  output  (correlation  function)  runs  from  0 
to  2 (K-l)  with  2K-1  non-zero  elements. 
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APPENDIX  B 


SIGNAL-TO-NOISE  RELATIONSHIPS  FOR  HF  MODEMS 


For  completeness,  we  include  this  appendix  which  relates  the 
various  signal-to-noise  ratios  used  in  this  report.  A 39- tone 
modem  is  specified  with  a frame  rate  of  22.5  ms,  and  an  integra- 
tion time  of  17.78  ms,  i.e.,  a guard  time  of  4.72  ms,  is  used. 

A complex  signaling  rate  of  3600  Hz  is  consistent  with  a digital 
implementation  using  a 64-point  FFT  for  the  above  tone  specifica- 
tions. The  modulation  format  for  the  preamble  design  is  based  on 
using  only  four  parallel  tones  followed  by  a wideband  pulse  sequence. 
Thus,  any  definition  of  SNR  must  also  be  suitable  for  preamble 
evaluation. 

In  all  cases  it  is  meaningful  to  obtain  results  in  terms  of 
the  ratio  of  the  total  received  power  to  noise  power  density.  This 
signal-to-noise  density  ratio  ( SNDR)  is  given  by  the  ratio  of  the 
total  received  power,  P,  to  the  noise  power  per  Hz  of  bandwidth, 
denoted  by  Nq,  as 


SNDR  = ~ (1) 

For  preamble  evaluation  in  Sections  2 and  3,  results  are  given  in 
terms  of  this  SNDR.  The  above  is  called  a signal-to-noise  density 
ratio  since  the  denominator  is  the  noise  density,  rather  than  the 
more  conventional  noise  power,  in  the  bandwidth  of  the  signal.  To 
avoid  the  ambiguity  which  results  from  different  bandwidth  defini- 
tions, one  may  define  the  signal-to-noise  ratio  per  information  bit 


(2) 
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where  R is  the  actual  data  rate  of  the  modem.  For  speech  trans- 
mission R “ 2400  b/s  and  thus  Eqs . (1)  and  (2)  can  be  related  by 

:£■  - [~  + 33.8]dB  - for  R - 2400  b/s  (3) 

N0  \ N0  / 

For  digitized  speech  the  modulated  signal  is  approximately 
contained  in  a bandwidth  of  2400  Hz,  and  thus,  E^/Nq  can  also  be 
viewed  as  the  ratio  of signal-to-noise  power  measured  in  a band- 
width of  2400  Hz. 

In  addition  to  obtaining  results  in  terms  of  P/Nq  and 
E,/Nq , which  are  easily  related,  it  is  customary  for  parallel  tone 
evaluations  to  present  results  in  terms  of  the  average  signal-to- 
noise  ratio  per  tone  defined  by 


ET  „ Etr 2] 

N0  E[n*]  + E[n?] 
c s 


(A) 


where  r is  the  tone  amplitude  and  n and  n are  the  cosine  and 

c s 

sine  noise  components  at  the  output  of  the  integrate-and-dump 
circuits.  This  SNR  is  also  closely  related  to  the  signal-to- 
noise  power  measured  in  a bandwidth  equal  to  that  of  the  modulated 
signal . 

The  components  of  Eq.  (4)  are  obtained  by  first  noting  that 
a multiple  tone  received  signal  can  be  written  as 

N 

y(t)  = Jl  £ r . cos(cu,  t + 0.)  + n(t)  (5) 

i=l  1 
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where  n(t)  is  assumed  to  be  a white  Gaussian  nois  >rocess 


normalized  such  that  E[n(t)n(t  + ->  ) ] = (Nq/2)6(t). 


The  cosine  and  sine  components  per  tone  are  given  by 


•ft 


yc  mTl  y ^ COS  “i*  dt  = ri  COS  6i  + n< 


(6) 


ft 


y - ~zr  P y(t)  sin  co.t  dt  = r,  sin  0.  + n 

s . i J 1 1 li 


Note,  the  multiple  tones  are  orthogonal  over  the  integration  time 
T,  which  is  less  than  a baud  interval  by  the  allowable  guard 


time. 


The  average  signal  energy  per  tone  is  given  by 


Et  - E[r?  cos2  0^  + r2  sin2  0^  * E[r2] 


(7) 


and  the  energy  per  noise  component  is  given  by 


E[n2  ] = E[n2  ] = N /2 
c . s . u 

i i 


(8) 


The  total  received  power  can  be  written  as 
NEm 


P = 


(9) 


where  N is  the  total  number  of  tones  and  T is  the  integration  time 
per  tone.  Thus,  we  can  relate  (1)  and  (2)  to  E^/Nq  by  the  following: 
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No  T \N0 


(10) 


N0  TR  ^N0 


(11) 


For  a conventional  16  tone  modem*  operating  at  2400  b/s,  and 
including  a Doppler  correction  tone  of  twice  the  amplitude, 

^effective  = » e9uat;i-on  (11)  gives  the  rela  konship  between 

Efe/NQ  and  Et/Nq  as 


2400 


(12) 


Eb  Et 

^ - .38  dB  - 16  tone  modem 

INq 


(13) 


For  an  ANDVT  modem  operating  at.  2400  b/s  with  two  Golay 
codes,  using  39  tones  and  an  integration  of  17.78  ms,  we  have 


N0  17.78  x 10"3x  2400  \N0 


(14) 


D.  Chase,  A Combined  Coding  and  Modulation  Approach  for  Com* 
munication  Over  Dispersive  Channels,"  IEEE  Trans . on 
Communications . Vol.  COM-21,  No.  3,  March  1973,  page  159. 
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or 


dB  - 39  tone  ANDVT  modem 


(15) 


Thus,  in  this  case,  a comparison  based  on  Ej,/Nq  is  essentially 
equivalent  (adjustment  is  only  .01  dB)  to  the  fair  comparison 
based  on  E^/Nq. 

If  the  ANDVT  modem  is  used  without  coding  a data  rate  of 
2400  bps  can  be  achieved  with  27  tones  which  are  keyed  every 
22.5  ms  with  one  of  four  phases  (2  bits/tone).  In  this  case 
the  relationship  between  E^/Nq  and  Ej/Nq  is  given  by 


27 

17.78  x l(f 3 x 2400 


(16) 


S _ 

N0 


1.99 


dB  - 27  tone  ANDVT  modem 


(17) 


Note,  by  comparing  Eqs.  (13)  and  (17),  without  coding  the  ANDVT 
modem  offers  a 1.61  dB  advantage  in  E^/Nq  when  compared  to  the 
16-tone  modem.  This  gain  is  due  to  the  fact  that  the  loss  due 
to  guard  time  is  smaller  for  the  ANDVT  modem  and  also  due  ;:o  the 
Doppler  tracking  algorithm  for  the  ANDVT  modem  which  does  not 
require  an  external  tone. 

The  ANDVT  speech  modem  results  are  usually  presented  in  terms 
of  E^/Nq  and  the  preamble  results  are  presented  in  terms  of  P/Nq. 
Equation  (3)  results  relate  these  two  measures.  The  results  for 
the  KG  synchronization  sequence  are  presented  in  terms  of  E^/Nq. 


B-5 


which  has  the  desirable  property  that  guard  time  and  code  redun- 
dancy does  affect  this  SNR  measure.  As  a consequence  of  the  low 
data  rate  used  for  the  KG  sync  sequence,  E^/Nq  is  considerably 
greater  than  E^,/Nq  for  this  case,  and  thus,  it  is  more  meaningful 
to  make  comparisons  based  on  P/Ng.  Equation  (10)  relates  P/Ng 
to  Etj,/Nq,  which  for  a 39-tone  preamble  and  integration  of  17.78  ms 
yields 


“ [E^,/Nq  + 33.4]  dB  - 39-tone  KG  preamble  (18) 

In  summary,  the  preamble,  KG  sync  sequence,  and  multiple  tone 
modems  can  be  compared  in  terms  of  P/Nq  but,  by  convention,  E^/Nq 
and  E^,/Nq  are  sometimes  used.  Equation  (3)  provides  the  appropriate 
conversion  for  the  ANDVT  speech  modem,  and  Eq.  (18)  provides  the 
conversion  for  the  KG  sync  results  given  in  Section  4. 
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